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ABSTRACT 

We investigate a spatially-flat cold dark matter model (with the matter density param- 
eter Q m = 0.3) with a primordial feature in the initial power spectrum. We assume 
that there is a bump in the power spectrum of density fluctuations at wavelengths 
A ~ 30 — 60/i _1 Mpc, which correspond to the scale of superclusters of galaxies. There 
are indications for such a feature in the power spectra derived from redshift surveys 
and also in the power spectra derived from peculiar velocities of galaxies. We study 
the mass function of clusters of galaxies, the power spectrum of the CMB temperature 
fluctuations, the rms bulk velocity and the rms peculiar velocity of clusters of galaxies. 
The baryon density is assumed to be consistent with the BBN value. We show that 
with an appropriately chosen feature in the power spectrum of density fluctuations at 
the scale of superclusters, the mass function of clusters, the CMB power spectrum, 
the rms bulk velocity and the rms peculiar velocity of clusters are in good agreement 
with the observed data. 

Key words: galaxies: clusters: general - cosmology: theory - dark matter - large- 
scale structure of Universe - cosmic microwave background. 
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1 INTRODUCTION 

The spatially-flat cold dark matter (CDM) model (matter 
density parameter Q m = 0.3, flatness being restored by a 
contribution from a cosmological constant Qa = 0.7) with a 
scale-invariant initial conditions has been a standard model 
in cosmology over last five years (see e.g. Ostriker & Stein- 
hardt 1995). This model successfully explains many large- 
and small-scale structure observations including the mass 
function and the peculiar velocities of clusters of galaxies. 
The flat cosmological model with the matter density pa- 
rameter Q m = 0.3 is also consistent with the observations 
of Type la supernovae at redshift z ~ 1 (Perlmutter at al. 
1999; Riess et al. 1998). 

In this paper we investigate the CDM model with a 
primordial feature in the initial power spectrum of density 
fluctuations. Adams, Ross & Sarkar (1997) have noted that 
according to our current understanding of the unification 
of fundamental interactions, there should have been phase 
transitions associated with spontaneous symmetry breaking 
during the inflationary era. This may have resulted in the 
breaking of scale-invariance of the initial power spectrum. 
Chung et al. (2000) studied an alternative mechanism that 
can alter classical motion of the inflaton and produce fea- 
tures in the initial power spectrum. They showed that if 
the inflaton is coupled to a massive particle, resonant pro- 
duction of the particle during inflation modifies the evolu- 
tion of the inflaton, and may leave an imprint in the initial 
power spectrum. The spectral features in the initial spec- 



trum may also be generated if the inflaton evolves through 
a kink in its potential (Starobinsky 1992; Lesgourgues, Po- 
larski & Starobinsky 1998; Gramann & Hiitsi 2000). 

A number of different non-scale invariant initial con- 
ditions have been recently used to analyze the cosmic mi- 
crowave background (CMB) data (see e.g. Kanasawa et al. 
2000; Barriga et al. 2000; Atrio-Barandela et al. 2000; Grif- 
fiths, Silk & Zaroubi 2000; Hannestad, Hansen & Villante 
2000; Wang & Mathews 2000). Griffiths, Silk & Zaroubi 
(2000) and Hannestad, Hansen & Villante (2000) showed 
that the CMB data favour a bump-like feature in the power 
spectrum around a scale of k = 0.004/iMpc -1 . Barriga et 
al. (2000) studied the step-like spectral feature in the range 
k ~ (0.06— 0.6)/iMpc~ 1 and found that such a spectral break 
enables a good fit to both the APM and CMB data. Atrio- 
Barandela et al. (2000) investigated the temperature power 
spectrum in the CDM models, where the power spectrum of 
density fluctuations at z ~ 10 3 was in the form Pik) ~ fc -1,9 
at wavenumbers k > 0.05/iMpc -1 . This power spectrum of 
density fluctuations was derived by Einasto et al. (1999) by 
analyzing different observed power spectra of galaxies and 
clusters of galaxies. Atrio-Barandela et al. (2000) found that 
this form of the power spectrum of density fluctuations is 
consistent with the recent CMB data. However, in this pa- 
per we examine the mass function of clusters of galaxies in 
the same model and find that for f2 m = 0.3, the number 
density of clusters is significantly smaller than observed. 

Suhhonenko & Gramann (1999, hereafter SG) studied 
properties of clusters of galaxies in two cosmological models 
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which rely on the observed power spectra of the distribution 
of galaxies. In the first model (hereafter model 1), the power 
spectrum of density fluctuations at z ~ I0 3 was in the form 
P{k) ~ k~ 2 at wavenumbers k > 0.05/iMpc" 1 . In the second 
model (hereafter model 2), the power spectrum contained 
a feature (bump) at wavenumbers k ~ 0.1 — 0.2ftMpc _1 
(A ~ 30 — 60/i _1 Mpc), which correspond to the scale of 
superclusters (see e.g. Einasto et al. 1997). SG examined 
the mass function, peculiar velocities, the power spectrum 
and the correlation function of clusters in both models and 
found that in many aspects the power spectrum of density 
fluctuations in model 2 fits the observed data better than 
the simple power-law model 1. This study suggested that 
probably at wavenumbers k ~ 0.05 — 0.2MVIpc _ , the power 
spectrum of density fluctuations is not a featureless simple 
power law. 

Fig. 1 shows the power spectrum of density fluctua- 
tions in the CDM model examined in this paper (see eq. 
(5) below). We assume that there is a feature (bump) in the 
power spectrum at wavenumbers k = 0.1 — 0.2/i _1 Mpc. The 
power spectrum in the CDM model with a scale-invariant 
initial power spectrum is also plotted. It is assumed that 
the density parameter Q m — 0.3 and the normalized Hubble 
constant h = 0.65. For comparison, we show in Fig. 1 the ob- 
served power spectra derived from the distribution of galax- 
ies in the APM, Stromlo-APM and Durham/UKST surveys 
(Baugh & Efstathiou 1993; Tadros & Efstathiou 1996; Hoyle 
et al. 1999). For the Stromlo-APM and Durham/UKST sur- 
veys, we present estimates for the flux-limited sample with 
P(k) = 8000/j~ 3 Mpc 3 in the weighting function (see Tadros 
& Efstathiou 1996; Hoyle et al. 1999 for details). There are 
indications for a similar bump in the power spectrum derived 
from the Stromlo-APM survey. On the other hand, there is 
no similar feature in the APM and Durham/UKST power 
spectrum. Hoyle et al. (1999) analyzed the power spectrum 
for different volume-limited and flux-limited samples drawn 
from the Durham/UKST redshift survey. There are indica- 
tions for a similar bump in the power spectrum measured in 
a volume limited sample with z max = 0.08 (see Hoyle et al. 
(1999) for details). 

We also examined the power spectrum derived from the 
distribution of galaxies in the SSRS2+CfA2 redshift survey 
(da Costa et al. 1994). There is a feature in the SSRS2+CfA2 
power spectrum at wavenumbers k = 0.1 — 0.2/i _1 Mpc. Sil- 
berman et al. (2001) studied the power spectrum of peculiar 
velocities of galaxies and found an indication for a wiggle in 
the power spectrum: an excess near k ~ 0.05/i Mpc -1 and 
a deficiency at k ~ 0.1ft Mpc -1 . This wiggle in the power 
spectrum is similar to the spectral feature studied in this 
paper. Most recently, there are indications for such a fea- 
ture in the preliminary power spectrum derived from part 
of the 2dF redshift survey (Percival et al. 2001). 

SG studied the power spectrum of clusters using N- 
body simulations and showed that the power spectrum of 
clusters in model 2 is in good agreement with the ob- 
served power spectrum of the APM clusters determined by 
Tadros, Efstathiou & Dalton (1998). SG investigated also 
the relation between the power spectrum of clusters and the 
power spectrum of matter fluctuations and found that in 
this model the relation between the cluster power spectrum 
and matter power spectrum is not linear at wavenumbers 
k > O.lftMpc -1 (see also Gramann & Suhhonenko 1999). 
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Figure 1. The power spectrum of density fluctuations in the 
CDM model examined in this paper (solid line) and in the stan- 
dard CDM model with a scale-invariant initial power spectrum 
(dotted line). In the models studied, Q m = 0.3 and h = 0.65. 
Filled circles, open circles and crosses show the power spec- 
trum of the galaxy distribution in the APM, Stromlo-APM and 
Durham/UKST surveys, respectively. For clarity, error bars are 
only shown for the APM data. 



In this paper we investigate the mass function of clus- 
ters of galaxies and the temperature power spectrum in the 
model with a bump in the power spectrum of density fluc- 
tuations at the scale of superclusters of galaxies. We also 
study the rms bulk velocity of galaxies and the rms peculiar 
velocity of clusters of galaxies in this model. The results are 
compared with observations. We examine the flat cosmolog- 
ical model with the density parameter Q, m — 0.3, the baryon 
density Qth 2 — 0.019 and the normalized Hubble constant 
h — 0.65 and h — 0.70. These values are in agreement with 
measurements of the density parameter (e.g. Bahcall et al. 

1999) , with measurements of the baryon density from abun- 
dances of light elements (O'Meara et al. 2001; Tytler et al. 

2000) and with measurements of the Hubble constant us- 
ing various distance indicators (Freedman et al. 2000; see 
also Parodi et al. 2000). The Hubble constant is given as 
Ho — 100ft km s _1 Mpc -1 . To restore the spatial flatness 
in the low density model, we assume a contribution from a 
cosmological constant J7a = 0.7. 

To study the mass function of clusters we use the Press- 
Schechter (Press-Schechter 1974, hereafter PS) approxima- 
tion. The transfer functions T(k) and the temperature power 
spectra are calculated using the fast Boltzmann code CMB- 
FAST developed by Seljak & Zaldarriaga (1996). The code 
CMBFAST has been modified to incorporate a primordial 
feature in the initial power spectrum. We assume that the 
initial fluctuations are adiabatic and that the initial density 
fluctuation field is a Gaussian field. In this case, the power 
spectrum provides a complete statistical description of the 
field. 

This paper is organized as follows. In Section 2 we study 
the mass function of clusters of galaxies and temperature 
power spectrum in our model, and compare the results with 
observations. In Section 3 we examine peculiar velocities of 
galaxies and clusters of galaxies. Discussion and summary 
are presented in Section 4. 
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2 THE MASS FUNCTION OF CLUSTERS AND 
CMB ANISOTROPIES 

Let us first consider the CDM model, where the power spec- 
trum of density fluctuations at z ~ 10 3 is in the form 



P(k) = AkS(k)T 2 (k) 



f AkT 2 (k), 
\P(k Q )(k/k y 



if k < ko; 
if k > ko, 



where fco = 0.05ft Mpc . Here, the initial power spectrum 
is defined as Pi n (k) = AkS(k) and T(k) is the transfer func- 
tion, which describes the modification of the initial power 
spectrum during the era of radiation domination. The func- 
tion S(k) describes the deviation of the initial power spec- 
trum from a scale invariant form P(k) ~ k. The normal- 
ization constant A is determined by the large-scale CMB 
anisotropy. This form of the power spectrum at wavenum- 
bers k > 0.05ftMpc _1 was derived by Einasto et al. (1999) 
by analyzing different observed power spectra of galaxies 
and clusters of galaxies. We denote the CDM model, were 
the power spectrum is in the form (1), as model 1. 

Fig. 2a shows the power spectrum of density fluctu- 
ations in model 1 for h = 0.65 and h = 0.70. We also 
show the power spectrum in the CDM model with a scale- 
invariant spectrum (S(k) = 1). In comparison with the stan- 
dard model, the power spectrum of density fluctuations in 
model 1 is depressed at wavenumbers k > 0.05/iMpc -1 . A 
similar break in the power spectrum of density fluctuations 
was analyzed by Atrio-Barandela et al. (2000). 

To calculate the angular power spectrum of the CMB 
temperature fluctuations, we used the code CMBFAST, 
which was modified to incorporate the function S(k) in the 
initial power spectrum. We examined the models with no 
reionization (optical depth r = 0). Fig. 2b shows the CMB 
power spectrum, AT; 2 = 1(1 + 1)Ci/2tt, predicted in model 
1. Here Ci = (a 2 m ) and a; m are the coefficients of the spher- 
ical harmonic decomposition of the CMB temperature field: 
AT(6,ip) = ^ aimYim(9, ip). Fig. 2b shows also the tem- 
perature power spectrum predicted in the standard CDM 
model. We see that lowering the amplitude of the power 
spectrum at wavenumbers k > 0.05/i -1 Mpc, lowers also the 
amplitude of the temperature power spectrum at multipoles 
I > 400. The height of the first peak is different in the models 
in Fig. 2b due to the change in Hubble's constant. 

Fig. 2b shows also the CMB power spectrum derived 
from the Boomerang (Netterfield et al. 2001) and from the 
Maxima- 1 (Hanany et al. 2000) experiments. We see that 
the temperature power spectrum in model 1 is consistent 
with the observed temperature power spectrum. The ampli- 
tude of the second acoustic peak at / ~ 500 in the observed 
temperature power spectrum is smaller than that predicted 
in the standard CDM model with £l m — 0.3. However, in the 
analyzes we have assumed that the spectral index n = 1 and 
the optical depth r = 0. Also, both CMB datasets have a 
calibration uncertainty and a beam uncertainty that are not 
included in the errors plotted in Fig. 2b. Netterfield et al. 
(2001) analyzed the CMB power spectrum in the standard 
CDM model in more detail and showed that this model is 
consistent with the observed CMB data, once the parame- 
ters 1 and t^O, and beam and calibration uncertain- 
ties are taken into account. 

To study the mass function of clusters we use the Press- 
Schechter (1974, PS) approximation. The PS mass function 
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Figure 2. (a) The power spectrum of density fluctuations in 
model 1 for h = 0.65 (solid line) and for h = 0.70 (dot-dashed 
line). The dotted line and dashed line show the power spectrum 
in the standard CDM model for h = 0.65 and h = 0.70, respec- 
tively, (b) CMB power spectrum in the same models. The data arc 
from the Boomerang (circles) and Maxima-1 (triangles) experi- 
ments, (c) The cluster mass function in the same models. Squares 
and triangles show the cluster mass function derived by Bahcall 
& Ccn (1993) and Girardi ct al. (1998), respectively. The circle 
represents the result obtained by White et al. (1993) 



has been compared with N-body simulations (Efstathiou et 
al. 1988; White, Efstathiou & Frenk 1993; Lacey & Cole 
1994; Eke, Cole & Frenk 1996) and has been shown to pro- 
vide an accurate description of the abundance of virialized 
cluster-size halos. In the PS approximation the number den- 
sity of clusters with the mass between M and M + dM is 
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given by 



n(M)dM 



2 p b 



da(M) 



7r M a 2 (M) dM 



5i 



° XP f2^(M)J dM ' 
(2) 

Here p b is the mean background density and St is the linear 
theory overdensity for a uniform spherical fluctuation which 
is now collapsing; St = 1.675 for £7o = 0.3 (Eke et al. 1996). 
The function a(M) is the rms linear density fluctuation at 
the mass scale M. We will use the top-hat window function 
to describe halos. For the top-hat window, the mass M is 
related to the window radius R as M = 47rp 6 P 3 /3. In this 
case, the number density of clusters of mass larger than M 
can be expressed as 



n c i(> M) 



n(M')dM' 



(2tt) 3 / 2 



da(r) 



2 (r) dr 



■ exp 



5 2 t 



2a 2 (r) 



J- (3) 



Fig. 2c shows the cluster mass function predicted in 
model 1 for ft — 0.65 and ft — 0.70. The mass function 
in the standard CDM model is also plotted. We investigated 
the cluster masses within a 1.5ft" 1 Mpc radius sphere around 
the cluster center. This mass Mi. 5 is related to the window 
radius R as 



R = 8.43a; 



Mi. 



6.99 x 10 14 no/i _1 M s 



(ft^Mpc). 



(4) 

Here the parameter a describes the cluster mass profile, 
M(r) ~ r", at radii r ~ 1.5ft -1 Mpc. Numerical simula- 
tions and observations of clusters indicate that the parame- 
ter a « 0.6 — 0.7 for most of clusters (e.g. Navarro, Frenk & 
White 1995; Carlberg, Yee & Ellingson 1997). In this paper 
we use the value a — 0.65. 

Fig. 2c also shows the mass function of clusters de- 
rived by Bahcall & Cen (1993, hereafter BC) and by Girardi 
et al. (1998, hereafter G98). BC used both optical and X- 
ray observed properties of clusters to determine the mass 
function of clusters. The function was extended towards the 
faint end using small groups of galaxies. G98 determined the 
mass function of clusters by using virial mass estimates for 
152 nearby Abell-ACO clusters including the ENACS data 
(Katgert et al. 1998). The mass function derived by G98 
is somewhat larger than the mass function derived by BC, 
the difference being larger at larger masses (see Fig. 2c). 
Reiprich & Bohringer (1999) determined the cluster mass 
function using X-ray flux-limited sample from ROSAT All- 
Sky survey. They determined the masses for different outer 
radii of the clusters and for a radius r — 1.5ft -1 Mpc their 
mass function agrees with that determined by BC. 

Let us consider the amplitude of the mass function of 
galaxy clusters at Mi. 5 = 4 ■ lO 14 ft _1 M0. For this mass, 
the cluster abundances derived by BC and G98 are n(> 
M) = (2.0 ± 1.1) • 10" 6 ft 3 Mpc" 3 and n(> M) = (6.3 ± 1.2) ■ 
10" 6 ft 3 Mpc~ 3 , respectively. By analysing X-ray properties 
of clusters, White, Efstathiou & Frenk (1993) found that 
the number density of clusters with the mass Mi. 5 ^ 4.2 ■ 
lO 14 ft _1 M is n(> M) — 4 • 10" 6 ft 3 Mpc" 3 . 

Fig. 2c shows that the cluster mass function in the 



standard CDM model is in good agreement with the ob- 
served data. But the number density of clusters in model 1 
is substantially lower than observed; for the mass Mi. 5 = 
4 • 10 14 ft _1 M©, the cluster abundance n(> M) = 1.5 • 
10" 7 ft 3 Mpc" 3 and n(> M) = 4.1 • 10" 7 ft 3 Mpc" 3 for ft = 
0.65 and ft = 0.70, respectively. 

In model 1, the amplitude of the power spectrum of den- 
sity fluctuations is depressed with respect to the standard 
CDM model for k > 0.05ftMpc _1 . Lowering the amplitude 
of the power spectrum of density fluctuations at wavenum- 
bers k > 0.05ft Mpc -1 lowers also the CMB power spectrum 
at multipoles I > 400. However, lowering the amplitude of 
the power spectrum of density fluctuations lowers also the 
cluster mass function, and as a result in model 1 the num- 
ber density of clusters is smaller than observed. One possi- 
bility to get rid of the last effect is to consider a bump in 
the power spectrum of density fluctuations at wavenumbers 
k ~ 0.1 - 0.2ft _1 Mpc. The cluster mass function for masses 
M ~ 10 14 - lO 15 ft _1 M is sensitive to the amplitude of the 
power spectrum at wavenumbers k ~ 0.2ft _1 Mpc, while the 
temperature anisotropy at the second acoustic peak is sensi- 
tive to the amplitude of the power spectrum at wavenumbers 
k ~ 0.05 - 0.1ft _1 Mpc. 

Thus let us now study a CDM model, where the power 
spectrum of density fluctuations at z ~ 10 3 contains a spe- 
cific feature at wavenumbers k ~ 0.1 — 0.2ftMpc _1 (A ~ 
30 — 60ft -1 Mpc) which correspond to the scale of superclus- 
ters: 

AkT 2 (k), 



P(k)=AkS(k)T 2 (k) = { 2M (fc/feo) 
P(fci), 

BkT 2 (k), 



if k < ko; 
if kg < k < fci; 
if fci < k < ki\ 
if k > fo, 

(5) 

where fco = 0.05ftMpc 1 , the spectral index m = 
log[P(fci)/P(fco)]/log[/ci/fe ] and B = P(fci)/(fc 2 T 2 (fc 2 )). 
This form of the power spectrum contains three free pa- 
rameters, which describe the bump in the power spectrum: 
fci, &2 and P(fci). The parameter k\ determines the be- 
ginning of the bump, the parameter fc 2 - the end of the 
bump, and the parameter P(ki) - the amplitude of the 
power spectrum for the bump. In this paper we examine 
the models where ki = O.lftMpc -1 , fc 2 = 0.2ftMpc _1 and 
P(fci) = 2500 - 3500ft~ 3 Mpc 3 . We denote the CDM model, 
where the power spectrum is in the form (5), as model 2. 

Fig. 3a shows the power spectrum of density fluctu- 
ations in model 2 for different values of the parameter 
P(fci). In the models studied, P(fci) = 2500,3000 and 
3500ft~ 3 Mpc 3 . The normalized Hubble constant ft = 0.65. 
In Fig. 3b, we show the function S(k), which describes 
the deviation of the initial power spectrum from the scale- 
invariant form. The function S(k) = 1 at wavenumbers 
k < O.OSftMpc" 1 , reaches the minimum at fci = O.lftMpc -1 
and then increases up to the wavenumber fc 2 = 0.2ftMpc _1 . 
At the minimum, S(ki) = 0.48 and S{ki) = 0.67 for 
P(fci) = 2500ft" 3 Mpc 3 and P(fci) = 3500ft~ 3 Mpc 3 , re- 
spectively. For wavenumbers k > 0.2ftMpc~ , we find that 
S(k) = 1.3 and S(k) = 1.9, respectively. We investigated 
also the CDM model with ft = 0.70 assuming that P(fci) = 
2500,3000 and 3500ft" 3 Mpc 3 . In this model, S(ki) = 0.38, 
S(k 2 ) = 1.0 and S(ki) = 0.52, S(k 2 ) = 1.4 for the param- 
eter P(fci) = 2500ft" 3 Mpc 3 and P(fci) = 3500ft" 3 Mpc 3 , 
respectively. 
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Figure 3. (a) The power spectrum of density fluctuations in 
model 2 for P(fci) = 2500/i" 3 Mpc 3 (solid line), for P(jfci) = 
3000/i" 3 Mpc 3 (dot-dashed line) and for P(fci) = 3500/i- 3 Mpc 3 
(dashed line). The dotted line shows the power spectrum in the 
standard CDM model. (Here h = 0.65.) (b) The function S(k) in 
the same models as in the panel (a). 



Fig. 4 shows the variance, a 2 (R), in the standard CDM 
model and in model 2 for different values of the parameter 
P(fci). The variance is given as a function of top-hat window 
radius R. For the masses Mi. 5 = W 14 h~ 1 M Q and M1.5 = 
10 15 /i _1 M Q , the window radius R = 5.8ft _1 Mpc and R = 
15.3ft -1 Mpc, respectively (see eq. (4)). 

Fig. 5 shows the cluster mass function as predicted in 
model 2 for the parameter P(fci) = 2500, 3000 and 3500/i~ 3 
Mpc 3 . The mass function of clusters in the standard CDM 
model is also plotted. Fig. 5a shows the results for h = 0.65 
and Fig. 5b for h — 0.70. In model 2, the mass function is 
steeper than that in the standard CDM model. At the same 
value of P(fei), the cluster mass function for h = 0.65 and 
h = 0.70 is similar. For comparison, we also show in Fig.5 the 
observed mass function of clusters of galaxies derived by BC 
and G98. In the models studied, the cluster mass function 
is consistent with the observed data. If h = 0.65, then for 
the mass Mi. 5 = 4 ■ 10 14 /i _1 Mq, the cluster abundance n(> 
M) = 1.4 ■ 10" 6 fe 3 Mpc- 3 and n{> M) = 3.5 • 10" 6 /i 3 Mp C - 3 
for P(fci) = 2500ft" 3 Mpc 3 and P(fci) = 3500ft~ 3 Mpc 3 , 
respectively. (If h = 0.70, n(> M) = 1.5 ■ 10" 6 ft 3 Mpc" 3 and 
n(> M) = 3.7 • 10~ 6 ft 3 Mpc" 3 , respectively.) 

Fig. 6 demonstrates the angular power spectrum of the 
CMB temperature fluctuations as predicted in model 2. We 
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Figure 4. The variance u 2 (R) in the same models as in Fig. 3, 
shown with the same line types. 
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Figure 5. The mass function of clusters as predicted in model 
2. (a) h=0.65. (b) h=0.70. The lines are specified similarly to the 
Fig. 3. Squares and triangles show the mass function of galaxy 
clusters derived by Bahcall & Cen (1993) and Girardi et al. (1998), 
respectively. 



also show the temperature power spectrum in the stan- 
dard CDM model with a scale-invariant power spectrum. 
Fig. 6a demonstrates the results for h — 0.65 and Fig. 6b 
for h = 0.70. In model 2, the amplitude of the temperature 
power spectrum at multipoles I > 400 is smaller than that 
predicted in the standard CDM model. Fig. 6 also shows the 
CMB power spectrum derived from the Boomerang (Net- 
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Table 1. Peculiar velocities in the different models. 
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Figure 6. The power spectrum of the CMB temperature in 
the same models as in Fig. 3, shown with the same line types, 
(a) h=0.65. (b) h=0.70. The observational data are from the 
Boomerang (circles) and Maxima-1 (triangles) experiments. 

terfield et al. 2001) and from the Maxima-1 (Hanany et al. 
2000) experiments. In the models studied, the temperature 
power spectrum is consistent with the observed temperature 
power spectrum. 

Therefore, in the model with a bump in the power spec- 
trum of density fluctuations at the scale of superclusters, the 
mass function of clusters and the temperature power spec- 
trum are in good agreement with the observed data. 



3 PECULIAR VELOCITIES 

The observed rms peculiar velocity of galaxy clusters has 
been studied in several papers (e.g. Bahcall, Gramann & 
Cen 1994, Bahcall & Oh 1996, Borgani et al. 1997, Watkins 
1997, Dale et al. 1999). Watkins (1997) developed a like- 
lihood method for estimating the rms peculiar velocity of 
clusters from line-of-sight velocity measurements and their 
associated errors. This method was applied to two observed 
samples of cluster peculiar velocities: a sample known as the 
SCI sample (Giovanelli et al. 1997) and a subsample of the 
Mark III catalog (Willick et al. 1997). Watkins (1997) found 
that the rms one- dimensional cluster peculiar velocity is 
265^7g 6 km s _1 , which corresponds to the three-dimensional 
rms velocity 459±j| km s" 1 . Dale et al. (1999) obtained 
Tully-Fisher peculiar velocities for 52 Abell clusters dis- 
tributed throughout the sky between ~ 50 and 200/i~ 1 Mpc. 



a 


F(fel) 

(/t- 3 Mpc- 3 ) 


Vcl 

(km s- 1 ) 


T 7 

(km s _1 ) 


0.65 


2500 


557 


273 




3000 


578 


273 




3500 


598 


273 


0.70 


2500 


579 


284 




3000 


599 


284 




3500 


619 


284 



They found that the rms one-dimensional cluster peculiar 
velocity is 341 ± 93km s _1 , which corresponds to the three- 
dimensional rms velocity 591 ± 161 km s _1 . 

To investigate peculiar velocities of clusters in our mod- 
els, we use the linear theory predictions for peculiar veloci- 
ties of peaks in the Gaussian field. The linear rms velocity 
fluctuation on a given scale R at the present epoch can be 
expressed as 

a v {R) =H f{Sl )a- 1 {R), (6) 

where /(fio) ~ fio i s the linear velocity growth factor in 
the flat models and Oj is defined for any integer j by 



"•-hi 



P(k)W z (kR)k 23+2 dk 



(7) 



Bardeen et al. (1986) showed that the rms peculiar velocity 
at peaks of the smoothed density field differs systematically 
from a v (R), and can be expressed as 



(8) 



SG examined the linear theory predictions for the pe- 
culiar velocities of peaks and compared these to the pe- 
culiar velocities of clusters in N-body simulations. The N- 
body clusters were determined as peaks of the density field 
smoothed on the scale R ~ 1.5h Mpc. The numerical re- 
sults showed that the rms peculiar velocity of small clusters 
is similar to the linear theory expectations, while the rms pe- 
culiar velocity of rich clusters is higher than that predicted in 
the linear theory. The rms peculiar velocity of clusters with a 
mean cluster separation d c i = 30/i _1 Mpc was ~ 18 per cent 
higher than that predicted by the linear theory. We assume 
that the observed cluster samples studied by Watkins (1997) 
and Dale et al. (1999) correspond to the model clusters with 
a separation d c i ~ 30/2T 1 Mpc (n c i ~ 3.7 ■ 10 _5 /i 3 Mpc -3 ) 
and determine the rms peculiar velocity of the clusters, v c i, 
as 



v ci = 1-18 CTp(R), 



(9) 



where the radius R — 1.5/i _1 Mpc. 

In the standard CDM model with Q. m = 0.3, we found 
that Vci — 582 km s _1 and v c i = 635 km s _1 for h — 0.65 and 
h — 0.70, respectively. Table 1 lists the rms peculiar velocity 
of clusters, v c i, in our models. The rms peculiar velocity of 
clusters is ~ 555 — 620 km s _1 , which is consistent with the 
observed rms peculiar velocity of clusters derived by Watkins 
(1997) and Dale et al. (1999). 

We also studied the rms bulk velocity for a radius r = 
60/i _1 Mpc, V 60 . The rms bulk velocity was determined by 
using equation (6). In the standard Q. m = 0.3 model, Vqo = 
273 km s _1 and Veo = 285 km s _1 for h = 0.65 and h = 
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0.70, respectively. Table 1 lists the rms bulk velocity, Veo, 
in our models. The bulk velocity is similar to that in the 
standard model. The observed bulk velocities are determined 
in a sphere centred on the Local Group and represent a 
single measurement of the bulk flow on large scales. The 
observed bulk velocity derived from the Mark III catalogue 
of peculiar velocities for r = 60/i~ 1 Mpc is 370 ±110 km 
s" 1 (Kolatt & Dekel 1997). Giovanelli et al. (1998) studied 
the bulk velocity in the SCI sample and estimated that the 
bulk flow of a sphere of radius r = 60ft~ 1 Mpc is between 
140 and 320 km s~ . In the models studied, the rms bulk 
velocity is ~ 275 — 285 km s _1 , which is consistent with the 
observed data. 



4 DISCUSSION AND SUMMARY 

In this paper we have examined a CDM model, where the 
power spectrum contains a specific feature (bump) at the 
wavenumbers k ~ 0.1— 0.2/i _1 Mpc, which correspond to the 
scale of superclusters of galaxies. We studied a flat cosmolog- 
ical model with the density parameter Q m = 0.3 and the nor- 
malized Hubble constant h — 0.65 and h = 0.70. The baryon 
density was assumed to be consistent with the standard big- 
bang nucleosynthesis (BBN) value. We investigated the mass 
function of clusters and the angular power spectrum of the 
CMB temperature fluctuations, assuming different values of 
the spectral parameter -P(fci) that determines the amplitude 
of the power spectrum for the bump. We found that the clus- 
ter mass function and the CMB power spectrum are in good 
agreement with the observed data if the spectral parameter 
P(fei) is in the range P(fci) = 2500 -3500ft -3 Mpc 3 . We also 
investigated the rms peculiar velocity of clusters and the rms 
bulk velocity for a radius r = 60/i _1 Mpc. In the models 
studied, the rms peculiar velocity of clusters is ~ 555 — 620 
km s _1 and the rms bulk velocity is ~ 275 — 285 km s _1 , 
which are consistent with the observed data. 

Therefore, in many aspects the CDM model, where the 
power spectrum contains a feature at the scale of superclus- 
ters of galaxies, fits the observed data. This model predicts 
that there is a bump in the correlation function of clusters 
at separations r ~ 20 — 35/i~ 1 Mpc (Suhhonenko & Gra- 
mann 1999). In the future, accurate measurements of the 
cluster correlation function at these distances can serve as a 
discriminating test for this model. 

What could be the origin of the primordial feature ex- 
amined in this paper (see Fig. 3b)? The standard inflation- 
ary prediction concerning the initial power spectrum of den- 
sity fluctuations is a simple power law: Pi n ~ k n . How- 
ever, for more than ten years, there has been some inter- 
est in models called broken-scale-invariant (BSI), predict- 
ing deviations from a power-law. These models generally 
involve, in addition to the usual inflation field, other (ef- 
fective) fields, driving successive stages of inflation or just 
triggering a phase transition (e.g. Starobinsky 1985; Kof- 
man, Linde & Starobinsky 1985; Kofman & Pogosyan 1988; 
Gottlober, Miiller & Starobinsky 1991). Recently, there have 
been efforts to implement double/multiple inflation in realis- 
tic supersymmetric contexts (Lesgourgues 1998). There has 
also been an attempt to generate spectral features via res- 
onant production of particles during inflation (Chung et al. 
2000). Clearly, further work is needed to explain the pri- 



mordial feature introduced in this paper in the context of 
different inflationary models. 
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